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Abstract 

In this work we report Monte Carlo simulations of a 2D Ising model, in which the statistics 
of the Metropolis algorithm is replaced by the nonextensive one. We compute the magnetization 
and show that phase transitions are present for q ^ 1. A q — phase diagram (critical temperature 
vs. the entropic parameter q) is built and exhibits some interesting features, such as phases which 
are governed by the value of the entropic index q. It is shown that such phases favors some 
energy levels of magnetization states. It is also showed that the contribution of the Tsallis cutoff 
is essential to the existence of phase transitions. 
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1 Introduction 



The nonextensive statistics is a generalization of the Boltzmann - Gibbs one and it is based on the 
nonadditive entropy [I] 

1 - T V 3 

S q = k -^p (q G ») (1) 

where g is the entropic index for a specific system, connected to its dynamics, as recently proposed 
[21 E]; Pi are probabilities satisfying YliPi = 1, A; is a constant, and lim g _>i S q = Sbg, where Sbg 
is the Boltzmann - Gibbs entropy. In this statistics, a system composed of two independent parts 
A and B, in the sense that the probabilities of the systems factorize, has the following pseudo- 
additivity (nonextensivity) property of the entropy [H [5] 

S q (A + B) = S q {A) + S q {B) + (l-q) S q (A) S q (B)/k. (2) 

This pseudo-additivity is related to the composability property of S q [6]. Since for any system S q > 
0, then q < 1 correspond to superadditivity (superextensivity) , q = 1 to additivity (extensivity) , 
and q > 1 to subadditivity (subextensivity) . Besides representing a generalization, S q , as much as 
Sbg, is positive, concave and Lesche-stable (Vg > 0). Recently, it has been shown that it is also 
extensive for some kinds of correlated systems in which scale invariance prevails [8] . 

In this paper we report some results of a Monte Carlo simulation of a 2D Ising model upon 
replacing the statistics of the Metropolis algorithm by the nonextensive statistics. From numerical 
calculations we compute the magnetization of the system, and built a q — phase diagram showing 
that, even for q ^ 1, exist phase transitions, in contrast to a previous work The text is organized 
as follows: In Section 2, we describe the equilibrium distribution of nonextensive statistics and the 
importance of the internal energy constrains. In Section 3, we describe the introduction of the 
nonextensive formalism into the Monte Carlo method. In Sections 4 and 5, we discuss the main 
results and describe the behavior of the critical temperature with the entropic index in a phase 
diagram (T c vs. q). 



2 Nonextensive statistics 

To calculate the equilibrium distribution, the above entropy, Eq.([T|), must be maximized pp. If the 
system is isolated, i.e., in a microcanonical ensemble, 

Eft = 1 ( 3 ) 

8=1 

the maximization yields equiprobability of states occupation. On the other hand, if the system is 
in contact with a thermal reservoir (canonical ensemble), it is necessary to add the internal energy 
constraints, which can be done according to three possible choices. The first one is [1] 

n 

Y,Pi £ i = U (4) 

i=i 

the standard definition of internal energy in which {e{\ are the eigenvalues of the Hamiltonian of 
the system. The second, as postulated in [ID] , is: 



2 



1=1 



(5) 



Both definitions presents some difficulties with the interpretation of some results [UJ. Thus, a third 
choice for the internal energy constraint was introduced as 
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2_,i=i 

which defined the escort probability (first introduced in Ref. |12j): 



P 



(?) 



p 



4=1 P. 



The maximization of S q in that case yields for the probability distribution 



Pi = -7T 



(l-q)P (gj - C/g) 

2^7 = 1 Pi 



V(l-9) 



where 
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is the nonextensive partition function and f3 a Lagrange multiplier. After some algebraic manipu 
lations [UJ it becomes: 
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Pi 
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(10) 



and 



zL 



where 



n 

j=i 



and is the generalized exponential, which has the following the property: 



[l_(l_g)^ e .]l/(l-9) 



e" x , if 1 - (1 - q)x > 0; 
0, if 1 - (1 - q)x < 0. 



(11) 



(12) 



(13) 



known as the Tsallis cutoff procedure. A detailed discussion about the role of constraints within 
the nonextensive statistics was done by Tsallis et al [IT], but recently it has been showed by Ferri 
et al [13] the equivalence of all these formulations of internal energy constraints. In spite of that, 
in this work, to avoid misunderstanding, we choose the normalized internal energy form. 

The thermal equilibrium in the nonextensive statistics is still an open issue due to the definition 
of the physical temperature [IH [151 (13 El 13 13 [20]. Thus, differently from some authors [9], in 
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our approach the parameter j3' is assumed to be the physical temperature, i.e., 0' = (kT) . The 
validity of this choice was first shown experimentally [21] . and latter theoretically [2"1 [3"1 1221 123j for 
manganites. 



3 Monte Carlo simulations of a 2D Ising model using nonextensive 
statistics 

In this section we are going to discuss the modification of the Metropolis method for the nonexten- 
sive statistics, considering a ferromagnetic 2D Ising with first-neighbors interaction. The Hamilto- 
nian is given by 



Ti 



-J s i s j 

(ij) 



(14) 



where (ij) denotes the sum over first neighbors on a square lattice of size N = L x L, Si = ±1 
and J > (ferromagnetic interaction). We proceed the single flip Monte Carlo calculations [24] to 
obtain the magnetization of the system, however we have changed the usual statistical weight to: 
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^before 



/kT 



(15) 



or, in other words, the ratio between the escort probabilities before and after the spin flip. Since 
this quantity is a ratio, the normalization factor of the escort probabilities, i.e., the generalized 
partition function, Eq. ljlip . cancels and the weight calculated can be written as the ratio between 
the generalized exponentials raised to the entropic parameter q. It is important to emphasize that 
w is the quantity that will be compared to a random number in the Metropolis algorithm (see 
appendix for details on the MC procedure used) . It is also important to note that the Tsallis cutoff 
procedure, Eq. (fT3j) . must be taken into account, i.e., it must be included into w to avoid complex 
probabilities. The simulations were done with the entropic parameter q £ [0, 1]. Lattice size were 
L = 8, 16, 24, 32, with periodical boundary condition imposing that sl+i = s±. 



4 Results and Discussion 

The most probable normalized magnetization, m = M/L 2 , was obtained after 5 x 10 5 Monte Carlo 
steps and are shown on Fig.([Tj) for q = 0.4, q = 0.7 (which are representative results for q < 0.5 and 
q > 0.5, respectively), and for q = 1. One can observe that for q > 0.5 there are strong influences 
of the lattice size on the shape of the magnetization curve and on the critical temperature, T c . In 
addition, the magnetization drops smoothly to zero close to T c due to the thermal fluctuations. On 
the contrary, for q < 0.5, there are no dependence of those quantities on the lattice size, and the 
magnetization changes suddenly at T c , from m = 1 to m = 0. In other words, there are no thermal 
fluctuations in this case and the magnetization works like a microcanonic two-level system. 

This behavior, for q < 0.5, is simple to be understood. At low temperatures (for instance 
T = 1.5 and q = 0.4), the first Monte Carlo steps lead the magnetization to m = 1, i.e., to the 
ground state, as expected (due to the low temperatures). Then the subsequent Monte Carlo steps 
attempt to invert the spin, but it fails because is energetically unfavorable. Then the Metropolis 
algorithm takes place; as e a ^ ter = 4 J and e ^ ore = —4 J, therefore 1 — (1 — q) j5' e a J ter < for 
all T < 4 J (1 — q). So, considering the cutoff, Eq, (fT3|) . p^ ter = and then w = 0. Since the 
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Figure 1: Normalized most probable magnetization vs. temperature. For q = 0.4 it can be seen 
that the magnetization has no dependence on the lattice size, dropping form m = 1 to m = 
suddenly at T c . It happens due to the contribution of the Tsallis cutoff. However, for q = 0.7 the 
magnetization depends on the lattice size and smoothly goes to zero close to T c as a second order 
phase transition. For q = 1 it recovers the well known result. 

Metropolis algorithm flips energetically unfavorable spins if the random number is smaller then w, 
for T < 4 J (1 — q) those spins never flips (w = 0), keeping the magnetization at m = 1; in other 
words, in the ground state. This situation persists up to T c = 4 J (1 — q), where the cutoff for r 
is no longer satisfied and the thermal fluctuation can therefore acts. However, the spins are already 
quite warm and the magnetization drops suddenly to zero, i.e., to a equiprobable state. A similar 
behavior was already found describing the generalized Brillouin function |22j. 

To determine the critical temperature we must take the thermodynamic limit (L — > oo). In 
Fig. (J2J) we plot the critical temperature as a function of the inverse lattice size for different values 
of the entropic parameter, taking the limit {L~ l — ► 0). Notice that the critical temperature for 
q = 1 tends to Onsager result and, as explained above, it does not depend on the lattice size for 
q < 0.5. Similar independency were found in different systems [25]. Also, the slope of the curves 
changes with the entropic parameter q, suggesting a dependence of the critical exponents on the 
entropic parameter (for studies about this connection see for example [261 E3 EE])- 

With those critical temperature values we build a phase diagram shown in Fig. ([3]), i.e., T c as a 
function of q. It is quite interesting because, in contrast to previous works [9], we found that for 
the Monte Carlo simulations of a 2D Ising model in nonextensive statistics has phase transitions 
for q 7^ 1. It is clear that below the 4 J(l — q) line the system is in the ground state and then 
m = 1. Above this line there are two regions where the thermal fluctuation act: one above T c , 
i.e., in the paramagnetic regime and, consequently, in the equiprobable state; and the other regime 
lies between 4 J(l — q) line and T c when the magnetization assume values between and 1. It is 
interesting that the slope of the critical temperature, T c (q), does not changes abruptly with the 
increase of the entropic parameter. For q ~ 0.5 the slope changes smoothly indicating that the spin 
is not warm enough and pass to a thermal distribution region before the equiprobable state. 
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Figure 2: Taking the thermodynamical limit of the system, one can notice that the behavior is 
quite different for q < 0.5 and for q > 0.5. As can be seen, for q = 0.4 the critical temperature 
has no dependence on the lattice size. It happens due to the contribution of the Tsallis cutoff. 
However, for q = 0.7 the lattice size dependence appears. For q = 1 the thermodynamical limit 
recovers the critical temperature of Onsager result (T c = 2 [arctan(l/\/2)] _1 ). 



5 Conclusions 

In this work, we studied a ferromagnetic 2D Ising model with first-neighbor interactions through 
a Monte Carlo simulation in which the Metropolis algorithm was changed to the nonextensive 
statistics. Magnetization as a function of the temperature for different values of q were evaluated 
and we found phase transition for q ^ 1, in contrast to a previous work [9]. This results arises due 
to the definition of the physical temperature. In addition, we also have showed the contribution of 
the Tsallis cutoff is of great importance and rules the phase transition for q < 0.5. 

The authors acknowledge S.M.D. Queiros, C. Tsallis and R. Toral for their comments. We would 
like to thanks the Brazilian funding agencies CNPq and CAPES. DOSP would like to thanks the 
Brazilian funding agency CAPES for the financial support at Universidade de Aveiro at Portugal. 



Appendix 

Each Monte Carlo step can be resumed as the following 

1. Compute the interaction energy of a given spin i— th of the lattice with its neighbors £? e f ore _ 
Y^j=i £ ij- After that, change the state of this spin and compute again its interaction energy, 

after T r after , before , ,1 1 r , , 

£j . If e i < £ i , accept the change of state; 

2. If the energy is not lower, using Eq. (|13j) . compute w. Compare this quantity to a number 
that belongs to the interval [0, 1] generated randomly. Being this random number smaller 
then w then accept the change of state, otherwise not. 

As can be seen, this is the ordinary Metropolis algorithm in which the probability of state was 
changed from the Boltzmann weight to the Tsallis factor, Eq.([7]). 
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Enlropic parameter q 



Figure 3: Three different regions can be seen in this q — phase diagram. The first one is the 
ferromagnetic ground state (m = 1) which is kept from zero Kelvin until T c = 4 J(l — q), due to 
the contribution of the Tsallis cutoff. The second region has influence at the thermal distribution 
(0 < m < 1). The third is the equiprobable paramagnetic state (m = 0). The error is smaller then 
the the size of those dots. 
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